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ABSTRACT 

We investigate the stability of strange quark matter and the properties of the corresponding 
strange stars, within a wide range of quark mass scaling. The calculation shows that the re- 
sulting maximum mass always lies between l.SAf© and LSAf© for all the scalings chosen 
here. Strange star sequences with a linear scaling would support less gravitational mass, and 
a change (increase or decrease) of the scaling around the linear scaling would lead to a higher 
maximum mass. Radii invariably decrease with the mass scaling. Then the larger the scaling, 
the faster the star might spin. In addition, the variation of the scaling would cause an order 
of magnitude change of the strong electric field on quark surface, which is essential to sup- 
port possible crusts of strange stars against gravity and may then have some astrophysical 
implications. 
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1 INTRODUCTION 

In studying the equation of state (EOS) of ordinary quark matter, 
the cruial point is to treat quark confinement in a proper way. Ex- 
cept the conventional bag mechanism (where quarks are asymptot- 
ically free within a large bag), an alternative way to obtain con- 
finement is based on the density dependence of quark masses, then 
the proper variation of quark masses with density would mimic the 
strong interaction between quarks, which is the basic idea of the 
quark mass-density-dependent model. 

Originally, the interaction part of the quark masses was as- 
sumed to be inversely proportional to the density (Fowler et al. 
1981; Chakrabarty 1991; Chakrabarty 1993; Chakrabarty 1996), 
and this linear scaling has been extensively applied to study the 
properties of strange quark matter (SQM). However, this class of 
scaling is often criticized for its absence of a convincing deriva- 
tion (Peng 2000). Then a cubic scaling was derived based on the 
in-medium chiral condensates and linear confinement (Peng 2000). 
and has been widely used afterwards (Lugones & Horvath 2003; 
Zheng et al. 2004; Peng et al. 2006; Wen et al. 2007; Peng et al. 

2008) . But this deriving procedure is still not well justified since 
it took only the first order approximation of the chiral conden- 
sates in medium. Incorporating of higher orders of the approxima- 
tion would nontrivially complicate the quark mass formulas (Peng 

2009) . In fact, there are also other mass scalings in the litera- 
tures (Dey et al. 1998; Wang 2000; Zhang et al. 2001; Zhang & 
Su 2002; Zhang & Su 2003). 

Despite the big uncertainty of the quark mass formulas, this 
model, after all, is no doubt only a crude approximation to QCD. 
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For example, the model may not account for quark system where 
realistic quark vector interaction is non-ignorable. However, we can 
not get a general idea of how the strong interaction acts from the 
fundamental theory of strong interactions in hand, i.e. QCD. Until 
this stimulating controversy is solved, we feel safe to take the prag- 
matic point of view of using the model. This work does not claim 
to answer how Nature works. However, it may shed some light on 
what may happen in interesting physical situations. In this respect, 
the quark mass-density-dependent model has been, and still is, an 
interesting laboratory. 

The aim of the present paper then, is to study in what ex- 
tent this scaling model is allowed to study the properties of SQM. 
To this end, we treat the quark mass scaling as a free parameter, 
to investigate the stability of SQM and the variation of the pre- 
dicted properties of the corresponding strange stars (SSs) within a 
wide scaling range. Furthermore, we try to demonstrate the general 
features of SSs related to astrophysics observations, whatever the 
value of the free parameters. 

The paper is organized as follows. In Section 2 we describe 
the formalism applied in calculating the EOS of the SQM in the 
quark mass-density-dependent model. In Section 3 we present the 
structure of the stars made of this matter, including mass-radius 
relation, spin frequency, electric properties of the quark surface. 
Finally in Section 4 we address our main conclusions. 



2 THE MODEL 

As usually done, we consider SQM as a mixture of interacting 
u, d, s quarks, and electrons, where the mass of the quarks niq 
(q = u, d, s) is parametrized with the baryon number density rib as 
follows: 
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m.qo + mi = mqo + 



C 



(1) 



where C is a parameter to be determined by stability arguments. 
The density-dependent mass includes two parts: one is the orig- 
inal mass or current mass mqo, the other is the interacting part mi. 
The exponent of density x, i.e. the quark mass scaling, is treated as 
a free parameter in this paper. 

Denoting the Fermi momentum in the phase space by i/i (i = 
It, d, s, e~), the particle number densities can then be expressed as 



rii = Qi 



and the corresponding energy density as 
Jo 



(2) 



(3) 



The relevant chemical potentials ^d, fJ's, and fie satisfy 
the weak-equilibrium condition (we assume that neutrinos leave the 
system freely): 



+ Me = Md, ^J■d = ^J■3■ 

For the quark flavor i we have 

de . dsi di^i 



(4) 
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We see clearly from Equ. JS} that since the quark masses are density 
dependent, the derivatives generate an additional term with respect 
to the free Fermi gas model. 
For electrons, we have 



The pressure is then given by 

p = — £ + 'y^^fiiUi 



(7) 



-fio + rii 



drrij ^ 
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dnb ^ \mj 
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(8) 



with 0,0 being the free-particle contribution 

^0 = 



487r2 



-|-3m^ arcsinh ( — - J 

\miJ _ 



(9) 



The baryon number density and the charge density can be 
given as: 

1 , 



rib = -{nu + rid + ris) 



Qq — Q riw 



1 1 



(10) 
(11) 
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Figure 1. The .stability window of the SQM at zero pressure with the quark 
mass scaling parameter x = 1/10,1/5,1/3,1,2. The stability region 
(shadow region), is where the energy per particle is lower than 930 MeV 
and two-flavor quark matter is unstable. 



The charge-neutrality condition requires Qq = 0. 
Solving Equs. (Hi, Go}, O, 

we can determine riu, rid, ris, 
and rie for a given total baryon number density rib. The other quan- 
tities are obtained straightforwardly. 

In the present model, the parameters are: the electron mass 
rrie = 0.511 MeV, the quark current masses rriuo, rritjo, rriso, the 
confinement parameter C and the quark mass scaling x. Although 
the light-quark masses are not without controversy and remain un- 
der active investigations, they are anyway very small, and so we 
simply take rriuo = 5 MeV, m^o = 10 MeV. The current mass of 
strange quarks is 95 ± 25 MeV accordi ng to the latest version of 
the Particle Data Group jYao et al. 2006h 

We now need to establish the conditions under which the SQM 
is the true strong interaction ground state. That is, we must re- 
quire, at P = 0,-B/A s; M(^^Fe)cV56 = 930 MeV for the 
SQM and EIA> 930 MeV for two-flavor quark matter (where 
Mi('^Y&) is the mass of ^''Fe) in order not to contradict standard 
nuclear physics. The EOS will describe stable SQM only for a set of 
values of (C, niao) satisfying these two conditions, which is given 
in Fig. 1 as the "stability window". Only if the (C, mso) pair is 
in the shadow region, SQM can be absolutely stable, therefore the 
range of C values is very narrow for a chosen mso value. As shown 
in Fig. 1, the allowed region decreases for decreasing value of x. 
When X = 1/10 it approaches to a very narrow area around C = 
199.1 MeVfm-3^. 

We then illustrate in Fig. 2 the density dependence of mi with 
the quark mass scaling a; = 1/10, 1/3, 1, 3. The calculation is done 
with mso = 95 MeV and C values corresponding to the upper 
boundaries defined in Fig. 1 (the same hereafter), that is, the sys- 
tem always lies in the same binding state (for each x), i.e, E/A = 
930 MeV. We presented those C values in the last row of the Ta- 
ble. 1 . Clearly the quark mass varies in a very large range from very 
high density region (asymptotic freedom regime) to lower densities, 
where confinement (hadrons formation ) takes place. It i s compared 
with Dey et al.'s scaling (dash-dotted) jDev et al. 19981) . 



3 RESULTS AND DISCUSSION 

The resulting EOSs of SQM are shown in Fig. 3 for all consid- 
ered models. Because the sound velocity v =| dP/dp \ should 
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Figure 2. The density dependence of mi with the quark mass scaling pa- 
ranieterx = 1/10, 1/3, 1, 3. The calculation is done with mso = 95 MeV 
and C values presented in the last row of the Table. 1 (see text for details). 
It is compared with Dey et al.'s scaling (dash-dotted line) (Dey et al. 1998). 



be smaller than c (velocity of light), unphysical region excluded 
by this condition has been displayed with scattered dots. For the 
X values chosen here, they have quite different behavior at low 
density, basically falling into two sequences. At small scalings 
(x — 1/10, 1/5, 1/3) the pressure increases rather slowly with 
density; while the curve turns to rapidly increase with density at 
relatively large x values (x — 1, 2, 3). They cross at e ~ 800 MeV 
fm then tend to be asymptotically linear relations at higher den- 
sities, and a larger x value leads to a stiffer EOS. Meanwhile we 
check also the stability of such a quark matter, since some EOSs in 
Fig 3 (a; = 1, 2, 3) are rather stiff for small pressures. We present 
in Fig. 4 the total pressure as a function of neutron chemical po- 
tential in quark matter for all considered models, and comparison 
with that of typical nuclear matter (obtained from the Brueckner- 
Hartree-Fock approach (Li et al. 2006)). We see clearly from the 
figure that the quark matter tends to be more stable than nuclear 
matter for all considered models. 

The behavior of EOSs would be mirrored at the prediction of 
mass-radius relations of the corresponding SSs, as is shown in the 
Fig. 5. For the first sequence, the maximum mass occurs at a low 
central density (as shown in Table. [TJ, so a higher maximum mass 
is obtained due to a stiffer EOS, and with the increase of x value, 
the maximum mass is reduced from USMq at a; = 1/10 down 
to 1.61 Mq at a; = 1/3; While we observe a slight increase of 
the maximum mass with x value for the second sequence: from 
1.56M0 at a; = 1 up to 1.62 AIq at a; = 3. Anyway the resulting 
maximum mass lies between 1.5Mq and I.8A/0 for a rather wide 
range of x value chosen here (0.1 - 3), which may be a pleasing 
feature of this model: well-controlled. To see the region of stellar 
parameters allowed by this model, we plot in Fig. 5 also the M(R) 
curves for lower boundaries defined in Fig. 1 for a; = 1/5, 1/3, 1 
with grey lines. 

The radii, on the other hand, decrease invariably with x 
value. Employing the empirical formula connecting the max- 
imum rotation frequency with the maximum mass and ra- 
dius of the static configuration jOourgoulhon et al. 1999I) . we 
get also the maximum rotational angular frequency S^max as 

rad s~^. As a result, a larger x value 
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Table 1. Calculated results for the gravitational masses, radii, central baryon 
densities (normalized to the saturation density of nuclear matter, no = 0.17 
fm^'^), and the maximum rotational frequencies for the maximum-mass 
stars of each strange star sequence. The calculation is done with rrigo = 95 
MeV and C values presented in the last row of this table. 
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Figure 3. The EOSs of SQM for all considered models. Unphysical region 
excluded by this condition has been displayed with scattered dots (see text 
for details). 



rotate at a frequency of /max = 2194 Hz. More detailed results can 
be found in Table. [T] 

In addition, the surface electric field could be very strong near 
the bare quark surface of a strange star because of the mass dif- 
ference of the strange quark and the up (or down) quark, which 
could play an important role in producing the thermal emission 
of bare strange stars by the Usov mechanism (Usov 1998; Usov 
2001). The strong electric field is also very crucial in forming a 
possible crust around a strange star, which has been investigated 
extensively by many authors (for a recent development, see Zdunik 
et al. 2001). Furthermore, it should be noted that this electric field 
may have some important implications on pulsar radio emission 
mechanisms (Xu et al. 2001). Therefore it is very worthwhile to 
explore how the mass scaling influences the surface electric field 
of the stars, and possible related astronomical observations in turn 
may drop a hint on what the proper mass scaling would be. 

Adopting a simple Thomas-Fermi model, one gets the Pois- 
son's equation (Alcock et al. 1986); 



2 > 0, 



(12) 



where z is the height above the quark surface, a is the fine-structure 
constant, and / {iiv^ h? (?) is the quark charge density inside 
the quark surface. Together with the physical boundary conditions 



results in a larger maximum spin frequency, SSs with a: = 3 can 
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Figure 4. The total pressure as a function of neutron chemical potential in 
SQM for all considered models, and comparison with that of typical nuclear 
matter (see text for details). 
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Figure 5. The mass-radius relations of SSs for all considered models. M(R) 
curves for lower boundaries defined in Fig. 1 with the quark mass scaling 
parameter x = 1/5, 1/3, 1 are presented with grey lines. Contours of the 
maximum rotation frequencies are given by the light grey curves (Gour- 
goulhon et al. 1999). 



0, dV/dz 0}, and the continuity of V al z — requires 
V{z = 0) = 3Vq/4, the solution for z > finally leads to 



V 



(forz > 0). 



(13) 



The electron charge density can be calculated as V'^ /{Siv'^h'^c'^) 
therefore the number density of the electrons is 



9V' 



n^{y/^V^z + 4r 
and the electric field above the quark surface is finally 



E = 



9V' 



(v^V',.z + 4)2 



(14) 



(15) 



which is directed outward. 

We see from Fig. 5 (take x = 1/3 for example) that although 
the electric field near the surface is about 10^* V cm~^, the out- 
ward electric field decreases very rapidly above the quark surface, 
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Figure 6. The electric field above the quark surface with the quark mass 
scaling parameter x = 1/3, 1, 2, 3. 



and at 2 ~ 10^ cm, the field gets down to ~ 5 x 10 V cm^ , 
which is of the order of the rotation-induced electric field for a typ- 
ical Goldreich-Julian magnetosphere. To change the mass scaling 
mainly has two effects: First, it affects a lot the surface electric 
field, and a small scaling parameter leads to an enhanced electric 
field. The weakening of electric field would be almost a order of 
magnitude large (from 10^^ V cm^^ to 10^** V cm^^), which may 
have some effect on astronomical observations. Second, a larger 
scaling would slow the decrease of the electric field above the quark 
surface. 



4 CONCLUSIONS 

In this paper, we investigate the stability of SQM within a wide 
scaling range, i.e. from 0.1 to 3. We study also the properties of 
the SSs made of the matter. The calculation shows that the result- 
ing maximum mass always lies between 1.5i\fQ and I.8A/0 for all 
the mass scalings chosen here. Strange star sequences with a linear 
scaling would support less gravitational mass, a change (increase or 
decrease) of the scaling parameter around the linear scaling would 
result in a higher maximum mass. Radii invariably decrease with 
the mass scaling. Then the larger the scaling, the faster the star ro- 
tates. In addition, the variation of the scaling may cause an order 
of magnitude change of the surface electric field, which may have 
some effect on astronomical observations. 



ACKNOWLEDGMENTS 

We would like to thank an anonymous referee for valuable com- 
ments and suggestions, and acknowledge Dr. Guang-Xiong Peng 
for beneficial discussions. This work was supported by the National 
Basic Research Program of China under grant 2009CB824800, 
the National Natural Science Foundation of China under grants 
10778611, 10833002, 10973002 and the Youth Innovation Foun- 
dation of Fujian Province under grant 2009J05013. 



REFERENCES 

Alcock, C, Farhi,E.. & Olinto, A. 1986 Astrophys. J. 310, 261 
Chakrabarty, S.; Raha, S., & Sinha, B. 1989 Phys. Lett. B, 229, 112 



Strange stars with different quark mass scalings 

Chakrabarty, S. 1991, Phys. Rev. D, 43, 627 
Chakrabaity, S. 1993, Phys. Rev. D, 48, 1409 
Chakrabarty, S. 1996, Phys. Rev. D, 54, 1306 

Dey, M., Bombaci, I., Dey, J., Ray, S., & Samanta, B. C. 1998, Phys. Lett. B, 

438, 123; erratum 1999, Phys. Lett. B, 467, 303 
Fowler, G. N., Raha, S. & Weiner, R. M. 1981, Z. Phys. C, 9, 271 
Gourgoulhon, E., Haensel, P., Livine, R., Paluch, E., Bonazzola, S., and 

Marck, J.-A., 1999, A&A 349, 851 
A. Li, G. F. Burgio, U. Lombardo, and W. Zuo, Phys. Rev. C74, 055801 

(2006). 

Lugones, G. & Horvath.J. E. 2003, Int. J. Mod. Phys. D, 12, 495 

Peng, G. X., Chiang, H. Q., Zou, B. S., Nmg, P Z., & Luo, S. J. 1999, Phys. 

Rev. C, 62, 025801 
Peng, G. X., Wen, X. J. & Chen, Y. D. 2006, Phys. Lett. B, 633, 313 
Peng, G. X., Li, A., & Lombardo U. 2008, Phys. Rev. C, 77, 065807 
Peng, G. X. 2009, private communication 
Usov, V. V. 1998, Phys. Rev. Lett., 80, 230 
Usov, V. V. 2001, ApJ, 550, L179 
Wang, P. 2000, Phys. Rev. C, 62, 015204 

Wen, X. J., Peng, G. X., & Chen, Y. D. 2007, J. Phys. G: Nucl. Part. Phys., 
34, 1697 

Xu, R. X., Zhang, B., & Qiao, G. J. 2001, Astroparticle Phys., 15, 101 
Yao, W.-M. et al. 2006, J. Phys. G: Nucl. Part. Phys, 33, 1 
Zdunik, J. L., Haensel, P, & Gourgoulhon, E. 2001, A&A, 372, 535 
Zhang, Y, Su, R. K., Ying, S. Q. Ying., & Wang, P. 2001, Europhys. Lett, 
56, 361 

Zhang, Y, Su, R. K. 2002, Phys. Rev. C, 65, 035202 
Zhang, Y, Su, R. K. 2003, Phys. Rev. C, 67, 015202 
Zheng, X. R, Liu, X. W., Kang, M., & Yang,S. H. 2004, Phys. Rev. C, 70, 
015803 



